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UNIVERSAL OPTIMALITY OF PATTERSON'S 
CROSSOVER DESIGNS^ 
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University of Waterloo, Indian Statistical Institute and Temple University 

We show that the balanced crossover designs given by Patterson 
[Biometrika 39 (1952) 32-48] are (a) universally optimal (UO) for the 
joint estimation of direct and residual effects when the competing 
class is the class of connected binary designs and (b) UO for the 
estimation of direct (residual) effects when the competing class of 
designs is the class of connected designs (which includes the connected 
binary designs) in which no treatment is given to the same subject 
in consecutive periods. In both results, the formulation of UO is as 
given by Shah and Sinha [Unpublished manuscript (2002)]. 

Further, we introduce a functional of practical interest, involv- 
ing both direct and residual effects, and establish (c) optimality of 
Patterson's designs with respect to this functional when the class of 
competing designs is as in (b) above. 

1. Introduction. Crossover designs (repeated measurement designs or 
change-over designs) in v treatments on n experimental units in p periods 
are useful in a broad spectrum of research areas, including agriculture [2], 
dairy husbandry [3] , bioassay procedures [4] , clinical trials [5] , psychological 
experiments [8] and weather modification experiments [17]. The advantages 
of the crossover design are its cost and the elimination of interunit variabil- 
ity. In the following, we assume that each treatment produces a direct effect 
in the period of its application and a residual effect in the subsequent period 
of its application. 

Williams [23] gave designs for p = v which were balanced in the sense 
that every paired difference of direct (residual) effects was estimated with the 
same precision. Patterson [18] gave combinatorial conditions for balance and 
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also gave a number of methods for construction of such designs when p <v 
and when n is as small as possible. Since p and n are small, these designs 
are very attractive to practitioners. All these designs had the property that 
no treatment immediately succeeds itself on the same subject. 

Hedayat and Afsarinijad [6] showed that when p = v, a balanced design 
is universally optimal (UO) (as defined in [9]) for estimation of the direct 
(residual) effects when the designs in the competing class are uniform on 
periods as well as subjects. Cheng and Wu [1] showed that these designs are 
UO for the estimation of residual effects when the competing designs may 
not be uniform over subjects or periods, but again no treatment succeeds 
itself on the same subject. Kunert [10] showed that when n = vt, a balanced 
uniform design is UO for direct effects if u > 3 and t = 1 or if u > 6 and 
t = 2. Hedayat and Yang [7] generalized this to the case where v >3 and 
t < {v — l)/2. The results of Kunert [10] and of Hedayat and Yang [7] were 
proved without any condition on the competing designs. However, there do 
not appear to be any available results on the optimality of balanced crossover 
designs when p <v. 

Cheng and Wu [1] also introduced what are called strongly balanced de- 
signs where each of the v"^ pairs of treatments occurs in consecutive periods 
for the same subject an equal number of times. They established some strong 
optimality properties for these designs. However, these designs require p = vt 
or vt-\-\ and also require n to be large. 

Kushner [14] gave a novel approximate design theory approach to obtain 
UO designs for arbitrary values of p and v. Further, Kushner [15] gave exact 
designs which are UO for direct effects for every pair (u,p) for some n. 

Kushner's results are very attractive because they do not put any condi- 
tions on the competing designs. Their main limitation is that the values of 
p or of n are large. Further, in almost all cases these optimal designs are 
nonbinary (on the subjects). An attractive property of the binary balanced 
designs is that they are optimal when the residual effects are negligible [9]. 

Some authors have obtained optimal designs under different models. Kunert 
and Martin [11] gave optimal designs under an interference model. Kunert 
and Martin [12] considered models with correlated errors. 

Kunert and Stufken [13] introduced a model where the residual effect of a 
treatment on itself is different from the residual effect when the treatment is 
followed by another treatment. An excellent review of the literature in this 
broad area up to 1996 is given by Stufken [22]. 

The balanced designs given by Williams [23] and by Patterson [18] are very 
attractive because they have a small number of periods and often involve 
a small to moderate number of subjects. These designs have been around 
for a long time and are generally believed to be efficient. However, precise 
optimality results are rather limited in nature. 



OPTIMALITY OF PATTERSON'S DESIGNS 



3 



In this paper we establish some strong optimahty properties of the Pat- 
terson designs. We first show that, within the class of binary designs, these 
designs are UO (in the sense of Shah and Sinha [21]) for the joint estimation 
of direct and residual effects. This is a very strong property because it im- 
plies UO for the estimation of the direct (residual) effects and a great deal 
more. For the rest of the paper we refer to UO as formulated in [21]. 

Next, we establish the UO property of these designs for the estimation 
of direct (residual) effects where the only restriction on a competing design 
is that no treatment immediately succeeds itself on the same subject. We 
also compute lower bounds for the efficiencies of these designs within the 
unrestricted class of competing designs and find that these are very high 
(0.99 or higher), giving rise to the speculation that, when the fully efficient 
designs as described by Kushner [15] do not exist, these designs might, in 
fact, be optimal for specific criteria, such as A-optimality. 

Further, we introduce an optimality function of practical interest and 
show that a Patterson crossover design is optimal for this functional, again 
with the restrictions on the competing designs that they are connected and 
that no treatment immediately succeeds itself on the same subject. 

In this paper we compare the information matrix for a Patterson design 
with the average (over permutations of treatment labels) of the information 
matrix for a competing design. This average has a form which is much sim- 
pler than the original matrix. Further, we first prove the optimality results 
for the model without period effects and then for the model with period 
effects. This works well because, when we introduce period effects, the infor- 
mation matrix for a Patterson design is unaltered, whereas, for a competing 
design, it is reduced by a nonnegative definite matrix. 

2. Preliminaries. Let us consider crossover designs where v treatments 
are arranged in p rows and n columns. The rows correspond to periods 
whereas the columns correspond to the subjects. A crossover design is said 
to be balanced if we have the following: 

(a) It is uniform over periods, that is, every treatment occurs t times in 
each period. 

(b) The design with subjects as blocks forms a balanced incomplete block 
design (BIBD). 

(c) The design with subjects as blocks and last period omitted also forms 
a BIBD. 

(d) Every ordered pair of distinct treatments occurs in consecutive pe- 
riods in units A times (any pair of identical treatments does not occur in 
consecutive periods). 

(e) In the set of t subjects receiving a particular treatment in the last 
period, every other treatment is applied A times in the first {p — 1) periods. 
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These conditions are equivalent to conditions I- VII given by Patterson. 
We shall call such designs Patterson designs. Existence of a Patterson design 
implies 

p<v, n = vt, X = t{p — l)/{v — 1), 

where t and A are positive integers. We shall assume that the parameter 
values p,v,n are such that a Patterson design exists. 

The crossover designs given by Williams [23] are balanced with p = v. 
However, one can often find balanced designs with p <v and where n is not 
too large. For v = A and p = 3, Patterson ([18], Figure 3) gives the following 
design with n = 12: 

abcbdadaccdb 
bcaabdacdbcd 
cabdabcdadbc 

Patterson [18] gave several methods of construction for balanced crossover 
designs with p <v. Our Table 1, above extracted from Table 1 in [18] giving 
designs available for small values of p, is of practical interest. 

Table 1 gives designs with the minimum value of n for given values of 
p and V when n < 60. Many more designs can be constructed using the 
various methods given by Patterson. Further, such designs are available with 
p = v = n when v is even and p = v = n/2 when v is odd [23]. Thus, this is a 
rich class of designs and it contains many designs of interest to practitioners. 

We exclude the case p = v = 2,n = 2t even though for these parameter 
values designs satisfying the combinatorial conditions exist. This is because 
in this case neither the direct nor the residual effects are estimable. 

Let d{i,j) be the treatment assigned to the jth subject in the ith period 
and let yij denote the response obtained from that subject in that period. 
We assume that the y^j's are uncorrelated with common variance o"^ and 

(2.1) E{yij) = + aj + 7j + Td{ij) + 

i = 1,2,.. . ,p;j = 1,2,. . . ,n;6d(o,j) = all j, where E{-) denotes the ex- 
pected value of the variable in the parentheses, fj, is the general mean, and 
a,7,r and 6 are the period, subject, direct and residual treatment effects, 
respectively. 

A crossover design is said to be connected if Tj — Tj/ and 6i — 6i' are es- 
timable for i^i! . All of Patterson's designs are connected. Let P denote the 
class of connected crossover designs using n (= vt) subjects for comparing 
V treatments in p periods, with the restriction that, in each column of a 
design, adjoining positions are occupied by distinct treatments. Further, let 
B denote the subclass of T) consisting of designs which are binary in the 
sense that no treatment is applied more than once to any subject. 

We define vectors Q:(p x l),7(?i x X),t{v x 1) and 8{v x 1) whose compo- 
nents represent the above effects: 
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p333344444555556666 
V 3 7 8 11 4 5 7 8 13 5 7 8 11 13 6 7 8 11 
n 6 21 56 55 4 20 14 56 52 10 21 56 55 39 6 42 56 22 



Let 



number of appearances of treatment i on subject u, 

number of appearances of treatment i on subject u 
in the first {p — 1) periods, 

number of appearances of treatment i preceded 
by treatment j on the same unit, 

number of appearances of treatment i in period k, 



0, 



i(fe-l) 



for k>2. 



We now define the frequency matrices N — (f^m), N — {niu), S — (•Sij), 
L = (^ifc), L = (lik)- Further, let diag(r) denote the v x v diagonal matrix 
whose elements are the replication numbers for the v treatments in the entire 
design. Also, let diag(f) denote the v x v diagonal matrix whose elements 
are the replication numbers for the treatments in the first (p — 1) periods 
only. 

The information matrix for {T,S,cx,'~f) is given by 



(2.2) 



T(t,5,q;,7) 



/diag(r) S 

diag(r) 

V ■ 



L 
L 

nln 



N 

pin/ 



where la denotes the identity matrix of order a and Jab denotes an a x 6 
matrix with all elements unity. (See equation (2.5) in [1].) 

The information matrix for {T,S,a) eliminating 7 is given by 



(2.3) T(T,<5,a|7) 



''diag(i 



V 



1, 



-NN* 



p 



-NN* 



diag(f ) 



-NN* 



P 







L- 


Inj 




1 


L- 


— NJnp 

P 


nlp 




--J 

Op 

p 



Here J^, is a /c x /c matrix with all elements unity. 
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When the period effects are ignored, the information matrix for (r, S) 
eliminating 7 is seen to be 



(2.4) 



X(r,5|7) 



diag(r) - -NN* 
P 



1. 



S - -NN* 

P 

-NN* 

P 



— I u 

Using ( w ) as a g-inverse of 

matrix for (t,(5) eUminating (a, 7) as 
I{T,S\a,f) 

(2.5) 



ml' 



diag(f) 

° ■^?"), we get the information 



/ 111* - — Nj 



■■I{t,S\j) 



n 



np 



LL* 



V 



n 
I. 

n 



-LL* 



np 
1 

np 



NJ„pL* 



NJ„pL* 



For a Patterson design, L = (tJ^jp), L = (0|tJt,^p_i), diag(r) =pilt,, diag(f ) = 
t{p - 1)1^. Further, NN* = p{t - \)\, +pU^, NN* = {{p-l)t-{p- 2) A)I^ + 
(p — 2)AJi, and S = A(Jt, — I^). We note that, for all designs in the design 
class 2?, the diagonal elements of S are all zeros. 

Without loss of generality, we arrange the subjects so that the first ni 
subjects have treatment 1 in the last period, the next n2 units have treatment 
2 in the last period and so on. 

This permits us to see that 













N* = N* - 





ln2 
















where 1/^ denotes an /i x 1 matrix with all elements unity. 
This gives us 



(2.6) 



NN* 
NN* 



NN* 
NN* 



and 

& - 0* + diag(ni,n2,...,n^), 



where = [0i , 02 j • • • > ^v] • Here, Oi is the sum of the nj columns of N cor- 
responding to the nj subjects where treatment i is in the last period. 
It is easy to verify that, for a Patterson design, 

= (t - A)!^, + AJ^,, NN*=p0, 
NN* = (p-l)0 and NN* = (p - 2)0 + 

For a design d, the information matrix for direct effects, residual ef- 
fects and period effects eliminating the subject effects, that is, Z{t,5,ol\^) 
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of (2.3), will be denoted by M^. Let d* denote a Patterson design. The 
matrix M^; for a Patterson design is given by 



p vt{p-l)^ vt{p - 1) ^ Q 



(2.7) Mrf. 



{v — 1) p{v — 1) 

^tip-^)„ t{p-l){pv-v^l) t{p-l) t 

7 \~ 7 ^ ' ^' ~ 

p(v — 1) p[v — 1) pv p 

p V 



Here H = 1^, — J^/v. 

Similarly, for del?, let denote the 2v x 2u information matrix for 
direct and residual effects eliminating the subject and period effects. We 
can write Cd as 

(2.8) IiT,S\c.,'y)=Cd= ^''''), 

where Cdu corresponds to the part for direct effects. Cdi2,Crf2i and 0^22 
are described similarly. Note that Cd2i = C^di2- Sometimes in the sequel we 
shall drop the suffix d. 

For the Patterson design d* , let = Cd*ij- Using (2.4), (2.5) and (2.8), 
we have 



(2.9) 



c;, = -^l<^^H and 

[v — 1) p[v — 1) 



t{p — l){pv — V — 1). 
p{v - 1) 



C22 — 7 7^ H. 



3. UO for joint estimation in the design class 13. In this section we shall 
show that d* is universally optimal (UO) for the joint estimation of direct 
and residual effects when the designs in the competing class are connected 
and are binary over subjects, that is, rnu = or 1. 

For formulations of UO one is referred to Kiefer [9] , Shah and Sinha [20] 
and Shah and Sinha [21]. Here, we shall use the formulation of Shah and 
Sinha [21], which may be described as follows. 

Let Cd denote a v x v direct effects information matrix (resp., v x v 
residual effects information matrix; or 2v x 2v joint direct-residual effects 
information matrix) of design d. Let g he a permutation of {1,2, ... ,v}, 
that is, g S S^, the symmetric group on {1,2,...,?;}. A design do having 
information matrix C^q is said to be UO in an appropriate design class if it 
minimizes every real valued function (j){C) (defined on the set of nonnegative 
definite matrices) that satisfies the following conditions: 

(1) (j){Cdg) = 0(Cd), where dg is the design obtained by permuting treat- 
ment labels according to g. 

(2) Cd >Cf^ 'A(Cd) < 0(Cj), where d and / are any two designs. 
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(3) 4>i^WgCdg) < 4>{^d)-, where Wg are all rational weights satisfying 
^Wg = l. Here g runs over all permutations in S^. 

It may be noted that every convex functional satisfies (3). This formula- 
tion of UO is an extension of Kiefer's original formulation [9], in the sense 
that the condition of convexity is replaced by the slightly weaker condition 
(3) above. See [20] and also [21] for a discussion of this. 

A sufficient condition due to Shah and Sinha [21] for dp to be UO is 

(3.1) ^WgCdg<Cd(, for every d 

In (3.1) the Wg^s can be any specific set of weights (which may depend 
upon d). In the sequel we will use (3.1) when Wg = l/v\ for all g € Sv 

Let Mrfg denote the matrix obtained from by permuting treatment 
labels according to g. We shall first show that 

(3.2) M.d*=Y,mdg/v\. 

9 

To show this, we shall state the following lemma which is easily established. 

Lemma 3.1. Let A be a k x k matrix and let g denote a permuta- 
tion on {1,2,..., A;} for which the permutation matrix is Pg. Then A = 
P^APg/A;! is a completely symmetric matrix with diagonal and off- 
diagonal elements a and b, respectively, given by 

a = '^^aii/k and b= — a^^ / k{k — \). 

Here s = J2i J2j (^ij sum of all elements of A. 

We now consider the various submatrices of for a binary design and 
show that, for each of these, the average over the permutations of treatment 
labels equals the corresponding expression for d* . 

For any binary design, the iih. diagonal element of NN* is X)i = 
^i^u'^iu = fii the replication number for the ith treatment. (This does 
not hold for a nonbinary design.) The average of over all permutations is 
pt, which is the replication number for d* . Further, for a binary design the 
iih. diagonal element of is rij. The average of ni over all permutations is 
t = n/v, the common diagonal element of for d* . 

Let 5 be a permutation on {1,2, . . . , ;;}. We note that the matrix M^^ is 

/Pg 

given by M.dg = Q^ M^Qg, where Q = o o 

\ I 

Using Lemma 3.1 and the expressions given in (2.3), one can easily verify 
that, if the design is binary and is connected for each of diag(r) — NN*/p, S — 
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NN*/p and for diag(f ) — NN*/p, the average over all g is indeed the corre- 
sponding expression for d* . We note here that, for the first two of these, the 
sum of all the elements is zero, whereas, for the third one, it is vt{p — l)/p- 

The assumption of connectedness is crucial here. If the design is not con- 
nected, NN* and/or NN* would have a block diagonal form where the 
off-diagonal submatrices consist of zeros. 

We also note that the average for L is tJyp and the average for N is 
{p/v)Jyn- Thus, the average for L — NJ„p/p is 0. 

Finally, we note that for each of L and N the average over the v\ per- 
mutations gives the corresponding expressions for a Patterson design. This 
completes the proof of the assertion. 

Now we note that the adjustment for the period effects a is equivalent 
to computing the Schur complement. Thus, the Schur complement of M^;* 
is Cd*, whereas the Schur complement of M^g is Cdg- Here C^* and Cdg 
refer to the 2v x 2v joint direct-residual effects information matrix. Since 
Mrf* = J2g^dg/vl and since the Schur complement is a concave function 
[19], we get 

Cd^ >J2Cdg/vl 
9 

Using the sufficient condition (3.1), with weights Wg = l/v\,g G , we see 
that d* is UO w.r.t. any design d£ B. 

As shown in [16] and in [21], UO for the joint estimation of two sets of 
parameters is a very strong property. In particular, it implies UO for the 
estimation of each set of parameters. It is shown in [21] that the converse 
is not true. A design can be UO for the estimation of the direct effects, as 
well as for the residual effects. However, it might fail to be UO for the joint 
estimation of the two. 

4. UO for direct (residual) effects in the design class X>. We now con- 
sider the case where the competing class of designs is V, a class of designs 
that contains all binary designs. Initially we shall assume that there are no 
period effects. We shall relax this assumption subsequently. 

The Cd matrix for this case is given by the submatrices of the matrix 
which correspond to the direct and the residual effects. These are the com- 
ponents of the information matrix Cd for the direct and the residual effects 
ignoring the period effects. We write the expressions for these dropping the 
suffix d, 

Cn = diag(r) - NN'/p, 
(4.1) Ci2 = S-NNVp, C2i = C*2, 

C22 = diag(f-) - NNVp- 



10 



K. R. SHAH, M. BOSE AND D. RAGHAVARAO 



We also write below the expressions for these submatrices for d* , the Pat- 
terson design, 

p* _p _ vt{p-l) 

(4.2) = Crf*i2 = - p^^_ = (ci;2)*) 

*-^22 — '-^d*22 — 7 TT H Jt;- 

p(f — Ij pv 

Let C denote the average of C over all permutations of treatment labels. 
To describe the structure of C, we introduce some notation. We define the 
following: 



(4.3) * " 

/ = ^ (sum of UiuS for subjects with treatment i in the last period). 

i 

Using the expressions for NN* and NN* in terms of NN* and given in 
(2.6) and using Lemma 3.1, one can verify that the matrices Cii,Ci2 and 
C22 have the following structure: 

Cii: Diagonal element is {p^vt — [5)/pv, 

off-diagonal element is —ijp'vt — f3)/pv{v — 1). 

C12: Diagonal element is —{(3 — l)/pv, 

off-diagonal element is (/3 — l)/pv{v — 1). 

C22: Diagonal element is {vt{p'^ — p — 1) — {f3 — sl))/pv, 

off-diagonal element is {(3 — 21 + pvt{2 — p))/pv{v — 1). 

To illustrate the nature of the computations, we consider C22- The sum 
of all the elements of diag(f) — NN*/p is vt{p — l)/p. The average of the 
diagonal elements of NN* = NN* - - 0* diag(ni, . . . , n^) is {(3-21 + 
vt)/v, whereas the average of the diagonal elements of diag(f) is t(p — 1). 
Use of Lemma 3.1 yields the expressions given above. 

From these we deduce that 

(A A) P{V - 1) P{V - 1) 

^ ' ^ pvt{p-l)-[P-2l)-t{v+p-l) ^ t{p-l) ^ 

p[v — 1) pv 

We shall now work with the information for direct (residual) effects ad- 
justed for residual (direct) effects and shall use condition (3.1) to show the 
UO property of Patterson designs. 
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We show below that both (C22) ^ and (C22) ^ exist. It is then easy to 
see that C11.22 = Cn — Ci2C^2^C2i for C* and for C are given by 

_vtip-l)( V 

(4.5) Cn.22 = ^^^^[p'tv-P 

ifi-l? 

[pvt{p - 1) - (/3 - 21) -t{v + p- 1)] 

We note that Cii,22 is the information matrix for the direct effects elim- 
inating the residual effects. Further, C11.22 is {J2g ^g/^^-) 11 ■22- 

We now show that each of C22 and C22 is nonsingular. We first consider 
C22- Since t{p — l)/pv 7^ 0, C22 is nonsingular iff the coefficient of H in the 
expression for C22 is nonzero. If this coefficient is zero, rank C22 is unity. 
We now recall that the {Si — (5j')'s are all estimable when the period effects 
are eliminated. These continue to be estimable when the period effects are 
ignored. Now, estimability of all {5i — 5j')'s implies that rank C22 >v — 1. 
This also implies that the sum of C22 over all permutations of treatment 
labels has rank at least {v — 1). We thus have rank C22 >v — 1 and, hence, 
the coefficient of H in the above expression for C22 must be nonzero. Prom 
the expression (4.2) for C22; it is clear that it is of full rank. 

We shall now show that C^j^ 22 ~ C11.22 is n.n.d. To see this, let kj denote 
the value of riiu in the jth of the rii subjects with treatment i in the last 
period. Similarly, let (3ij denote the contribution to /3 for that subject. We 
now note that (3 = ELiEjli Aj>^ = J2i=iJ2]Lihj- We note that / > vt. 
Further, for the j'th subject with treatment i in the last period, one riiu is 
lij and there are {v — 1) other nj^'s which add up to p — kj. Hence, 

Pij > lij + ip- kj) = kjikj -l)+p. 

It follows that 

(3ij - 2kj > ll - 3lij +p= {kj - l){lij -2)+p-2, 

Pij - kj > lij - 2lij +p= {kj -if+p-l. 

Since kj is 0, 1, 2 or greater than 2, we have 

(3ij - 2lij>p-2, (3ij -lij>p-l, (3ij>p. 

This gives f3 > pvt,(3 — I > vt{p — 1),(3 — 21 > vt{p — 2). Using the above 
relations, it is easy to see that C*;^ 22 — ^11.22, that is, C^;^ 22 ~ ^\\.22 is 
n.n.d. 

An application of (3.1) shows that d* is UO (compared with any design 
in V) if 

Cn.22 - ^{^g)\\.22lv\ is n.n.d. 
g 



|h and 
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Here (Cg)ii.22 is obtained by permuting treatment labels in C by and 
then computing the Schur complement (Cg)ii.22- In turns out that the result 
is the same as applying permutations g to the rows and columns C11.22, that 

is, (Cg)ii,22 = (Cii.22)g- 

Since the Schur complement C11.22 is a concave function [19], 

(EC.) >E(C,)i1.22. 

\ 9 / 11.22 9 

Since C11.22 = Eg Cg)ii.22, we have 

Cli,22 > Cn.22 = ( E C J /v\ > ^{Cg)u.22/Vl 
\ 9 I 11.22 9 

Thus, d* is UO in T) for the model without period effects. 

We shall now introduce period effects and consider estimation of the direct 
and the residual effects eliminating the period effects. When we adjust for 
periods, C gets reduced by an n.n.d. matrix and, hence, 

C (adjusted for periods) < C (ignoring periods). 

This implies ([19], Section 3.13) 

C (adjusted for periods) 11.22 < C (ignoring periods) 11.22- 

We now show that, for d*, 

C* (adjusted for periods) 11.22 = C* (ignoring periods) 11.22- 

Note that, for d*, the expressions for given by (2.9) and (4.2) differ 
only in C22- For these two cases, the Moore-Penrose inverses are = 
f(p-f}(p;-Li) H and = j(^;rfygS;rT)H + ^(^Jt,, respectively. Since 

C11.22 = Cii — C12C22C21, and since HJ„ = 0, the result follows. 
We now have 

C* (adjusted for periods)ii.22 = C* (ignoring periods)ii.22 

> ^E(Cg (ignoring periods))^! 22 

^' 9 

> — E^Cg (adjusted for periods));^^ 22- 
^' 9 

This establishes the UO property of d* (relative to any design in V) for 
the estimation of direct effects. 

An analogous argument also works for the estimation of residual effects. 
We outline the relevant important steps here. The information matrix for 
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residual effects eliminating direct effects is now given by C22.11 = C22 — 
C2iC^]^Ci2, where Cfi is the Moore-Penrose inverse of Cn. The expressions 
for C^2.ii and C22.11 = (Eg Cg/f )22.ii are given by 

r.* _ tip-l)ipv -V -1) ^ , t{p-l) ^ vt{p-[) 

^^22.11 — 7 T\ "^^ 2? TT ' 

p{v — 1) pv p-^iv — 1) 

pvt{p-l) - {p-2l)-t{v+p-l) 



^22.11 



p{v-l) 



H 



pv p{v — l)(jP'vt — (5) 

Use of 13 > pvt,/3 - I > vt{p - 1) and 13 - 21 > vt{p - 2) yields 
C22 11 — C22.11- As in the case of direct effects, one can now show the UO 
for the estimation of residual effects. 

It should be noted that d* is not UO in the class V for the joint estimation 
of direct and residual effects. This has been shown in [21]. 

It is not known if a Patterson design using the smallest number of subjects 
for given values of p and v is UO for the estimation of direct (residual) effects 
in the whole class of competing designs with fixed values of p, v and n. Hence, 
we obtain a lower bound for the efficiency factor for the estimation of direct 
effects along the lines of [15]. 

First, we note that, for an approximate optimal design, the information 
matrix C11.22 given by [15] is 

u — 1 L p pip — i-jv ) 

i;^(Cii.22) can serve as a lower bound to the (j) value of the information 
matrix for any competing design. Further, we note that both C11.22 and 
Cii,22 given by (4.5) are multiples of the matrix H. Hence, a lower bound 
to the efficiency of a Patterson design for the estimation of the direct effects 
is given by 

1-— .1/(1 + ' 



p{pv — V — 1) J / V pv{p — l) 
On simplification this reduces to 

e* = A/{A + v), 

where A = v'^{p — l)'^{pv{p — 1) — p — v). 

In Table 2 we give the values of e* for the 18 designs given in Table 1 of 
Section 2. 

As remarked earlier, these are very high. Thus, while our optimality re- 
sults are restricted to the design class V, we see that Patterson designs have 
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very high efficiencies even in the unrestricted class. As remarked in [15], 
when p <v, no design can be optimal in the unrestricted class for each of 
direct and residual effects. Patterson designs are UO for each of direct and 
residual effects within the class T> and are highly efficient even when we do 
not put any restrictions on the class of competing designs. 

5. A functional of interest. In this section we establish an interesting 
optimality property of d* when D is the class of competing designs. 

Sometimes, one may wish to look for the best combination of direct and 
residual effects. Thus, we would like to compare the yield when treatment 
j is followed by treatment i with the yield when treatment j' is followed by 
treatment i' . This means that we wish to estimate all functions of the form 



6j' as precisely as possible. This leads to the minimization of 



(5.1) 



A = ^ Var (fj + 5j - fi> - 6j'), 



where the summation is over such that 7^ j' ■ 

Since the variances of estimable parametric functions are functions of C, 
we may write A as ^(C). 



Let C^ 



n 



) denote the Moore-Penrose inverse of C. It is easy 
to see from (2.5) that each of fi, or A has zero row (column) sums. 
Expressing A as 

A = ^{Var(fi - f./ ) + Var((5j - Sy ) + 2 Cov(fi - f ./ , 6j - 6.>)} 

and using the above property of fi, and A, one can show that 



(5.2) A = 2v 



{(^-l) + (^-2)2}K^a,,+^A, -2(^-1)^/ 



a 



where fl = (aij), A = (Pij) and = (Oij). We note that (5.2) may not hold 
if row (column) sums of each ft, and A are not zero. 



Table 2 

Ejficiency lower bounds for Patterson 



designs for direct effects 



p 


V 


e* 


P 


V 


e* 


P 


V 


e* 


3 


3 


0.993103 


4 


7 


0.999783 


5 


11 


0.999972 


3 


7 


0.998885 


4 


8 


0.999835 


5 


13 


0.999980 


3 


8 


0.999156 


4 


13 


0.999939 


6 


6 


0.999960 


3 


11 


0.999563 


5 


5 


0.999853 


6 


7 


0.999971 


4 


4 


0.999306 


5 


7 


0.999930 


6 


8 


0.999978 


4 


5 


0.999565 


5 


8 


0.999947 


6 


11 


0.999988 
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One can see that A can be written as A = A{C) = {J2ieC^^^~^^)^'^ ^ where 
C denotes the set of coefficient vectors for all contrasts included in (5.1). 

We shall first show that A viewed as a functional on C and denoted by 
A{C) satisfies the three conditions on (piC) given in the UO formulation in 
Section 3. 

To see that A is invariant under a permutation of treatment labels, we 
shall show that a permutation g changes a coefhcient vector 1 into a vector 
gil) = 1', where 1' can be seen to satisfy all the necessary constraints. Let 
Pg be the v x v matrix for a permutation g £ S^. Let = (^^ p^). Then 

Cdg = F*C,F„ giving C+ = F*C+Fg. Thus, 1*C+ 1 = l'*C+l', where 1' = 
Fgl. As 1 varies over £, 1' also varies over C. Thus, condition (1) holds. 
Next, if Ci > Ca, < and, hence, 1*C^1 < 1*C^1 for ah 1 G C. Thus, 
condition (2) also holds. 

To show that condition (3) holds, we first note that C"*" is a convex func- 
tion of C, that is, {J2 '^g^g)^ ^ ^gC^, where the Wg's are rational weights 
satisfying ^Wg = l. Thus, we have 



1 



<E-.Ei*c+i..2 

9 1 

= WgA{Cg) = J2 WgA{C) = A{C). 

9 9 

This completes the verification. If we take Wg = l/v\, we get ^(C) < A{C). 
Thus, to show that ^(C,^*) < A(Crf), it is enough to show that A(Crf*) < 
A{Cd). 

Again, we initially assume that there are no period effects in the model. 
We first express C^* and Cd as 



(5.3) 



a*H b*H 

b*H c*H + eJy 
adH bdH 

bdil CrfH + eJy 



where a* ,b* ,c* ,e,ad,bd and Cd obtained from (4.2) and (4.4) are 

^* ^ vtjp - 1) ^ vt{p - 1) 

V — 1 ' p{v — 1) ' 

^* ^ t{p-l){pv - V -1) ^ ^ t{p - 1) 
p{v — 1) ' pv ^ 
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ad - 



Cd 



{13-1) 



p{v — 1) ' p{v — 1) ' 

pvt{p -l)-{j3- 21) -t{v + p-l) 



Cd- 



p{v - 1) 

For simplicity of notation, we shall drop the subscript d from dd,bd and 
Since C is an average of nonnegative definite (n.n.d.) matrices, it is n.n.d. 



Also, 0(^0^) =0 and 



Other eigenvalues of C are obtained as 



follows. Let uhe a. {v X {) vector satisfying Hu = u. It is easy to verify that 
(cm) eigenvector of C if q = {(c — a) it \J{a + c)^ — 4(ac — b"^)}/!. The 
corresponding eigenvalue is {(a + c) ± w (a + c)^ — 4(ac — 6^)1/2 . Since there 



are {v — 1) orthonormal choices for u, each of {(a + c) it y (a + c)^ — 4(ac — 6^) }/2 
is an eigenvalue of C of multiplicity v — 1. 

Since Tj — Tj' and 6j — 6j' are estimable in the model eliminating period ef- 
fects, they are also estimable in the model ignoring period effects. Thus, rank 
C>2{v-l). 

If A = ac — 6^ < 0, C has a negative eigenvalue. If A = 0, rank C = v, 



which is also a contradiction. Thus, A > 0. Similarly, A* 
Direct calculations yield 



6*2 > 0. 



(5.4) 



c*H/A* -6*H/A* 
-6*H/A* a*H/A* + J^./ev^ 

cH/A_ -6H/A 
-6H/A aii/A + J^/ev^ ^ 

Since the last v components in each vector 1 are v — 2 zeros, +1 and —1 
(in some order), we can, in computing X]ie£ 1*^)^1, ignore the term 
in C*"*" and C"*". The v x v submatrices of the remaining matrix have zero 
row and column sums. Hence (5.2) is applicable to this matrix. 

We can thus express A* = A{Cd*) and A = A{Cd) as 



A* = 2v{v-1) 



{(.-l) + (.-2)V" +^ 



A* 



+ 2{v-l) 



A* 



and 



A-- 



If we write (3 
c as 



2v{v - 1) 
21 = 



{{v-l) + iv-2) 



=}^ + 2(.-l) 



A ^ A_ 
vt{p — 2) + X and I = vt + y, we can express a, b and 



X + 2y 
p{v — 1) 



6 = 6* 



x + y 
p{v - 1)'- 



piv-1)' 
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It was noted in Section 4 that j3 — 2l> pvt{p — 2) and I > vt. Hence, we 
have X > and y > 0. 

We may write A* = A(0,0) and A = A{x,y). To show that A(0,0) < 
A{x,y), we shall proceed as follows: 

Using the expressions for a, 6 and c given above, it can be seen that 

A{x,y) = {(cii + ci2x + ci2.y) / {c2i + 022^; + C23y - y^/t{p - l))}o-^, 
where 

Cll = 2pv{v — l){2pv'^ — Qpv^ + 6pv — + 2v — 3), 
C12 = ci3 = -^pv{v - l){v^ -2v + 2)/t{p - 1), 
C21 =tv{p — l){p^v — pv — p — v), 
C22 = —{2pv + v — 1), C23 = —2{pv — 1). 
It is easy to verify that 

^ t(p — 1) , 2 / / -, M 

^ = ~27 T\2 + C22X + C23y - y t{p - 

— Ij^ 

Since A > 0, it follows that C21 + 0222^ + £23^ — y'^/t{p — 1) > 0. Similarly, 
A* > imphes C21 > 0. 

It follows that A{x,y) — ^4(0,0) is strictly positive iff (C21C12 — c\\C22)^ + 
(c2ici3 — ciiC23)y + ciiy"^ /t{p — 1) is strictly positive. For p>'i,v>p, the 
co-efficients of x^y and are all seen to be strictly positive and, hence, 
A{x, y) - A(0, 0) > if (x, y) / (0, 0). When p = 2, we must have y = 0. We 
shall comment on this case later in this section. 

We have thus shown that A{x, y) > A{0, 0) when period effects are ignored. 

When we take period effects into the model, C^* in (5.3) gets reduced by 
( ch) which has no effect on A. For C we argue as follows. We first note 
that 

C (adjusted for periods) < C (ignoring periods). 

Since C is obtained by averaging Cg over all permutations, it follows that 

C (adjusted for periods) < C (ignoring periods), 

C+ (adjusted for periods) > C+ (ignoring periods). 

Since ^(C) = ^jl*C+l • cr^, it follows that adjustment for periods cannot 
decrease the value of ^(C). 

We now summarize the situation as follows. Here, "adj" means adjusted 
for periods and "ign" means ignoring periods. We have seen that 

A(C<i*(adj)) = A(C<i*(ign)), 

^(Cd*(ign))<^(Cd(ign)), deV, 

A(Crf(ign))<A(Crf(adj)). 
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These imply 

A(Crf.(adj))<A(Cd(adj)), dev. 

This completes the proof of optimality of d* for the functional ^(C) in the 
design class V. 

In the definition of ^(C) we could also permit i = i' , as this would only 
add comparisons of the type 5j — 5ji. We have already seen that, for the 
estimation of residual effects, d* is UO in V. Similarly, we could also permit 
3=f- 

It should be noted that the above proof was needed only when d is non- 
binary with y > 0. When the design is binary, or nonbinary with y = 0, the 
result follows from the UO property of d* (Section 3 of this paper; [21]) for 
the joint estimation of direct and residual effects. 
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